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Abstract 

We show that the quasifission paths predicted by the one-body dissipa- 
tion dynamics, in the slowest phase of a binary reaction, follow a quasistatic 
path, which represents a sequence of states of thermal equilibrium at a fixed 
value of the deformation coordinate. This establishes the use of the statis- 
tical particle-evaporation model in the case of dynamical time-evolving sys- 
tems. Pre- and post-scission multiplicities of neutrons and total multiplicities 
of protons and a particles in fission reactions of 63 Cu+ 92 Mo, 60 Ni+ 100 Mo, 
63 Cu+ 100 Mo at 10 MeV/u and 20 Ne+ i44,i48,i54g m at 20 MeV/u are repro- 
duced reasonably well with statistical model calculations performed along 
dynamic trajectories whose slow stage (from the most compact configuration 
up to the point where the neck starts to develop) lasts some 35 x 10~ 21 s. 

PACS numbers: 25.70.Gh, 24.75. +i, 25.70. Jj 

1 Introduction 

In the last years many experimental efforts have been devoted to the study 
of heavy-ion fission at beam energies below 10-20 MeV/u accompanied with 
the emission of light particles. Initial experiments involved the observa- 
tion of fission fragments emitted in coincidence with neutrons |I|,@,§J or 
charged-particles 0|],§,[?|,|]|| . More recently, combined coincidence neutron 
and charged-particle data became available [^,|T!],|I2|,[nj . The most impor- 



tant result of these experiments consists in the estimation of fission times 
which appear to be 10 3 -10 4 times longer than the characteristic time of the 
single-particle motion. This fact implies that thermal equilibrium is being 
established over the intrinsic degrees of freedom at each value of the fission 
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coordinate Q and validates the use of the Rayleigh dissipative function tech- 



nique to treat the coupling of Q with the particle degrees of freedom |L4 . 

Heavy-ion fission reactions accompanied by particle emission are used as a 
testing ground of the dissipation mechanisms of large-scale collective motion 
in hot nuclei (for a review see Refs. [[i5i[l6| ). The main reason is that this 
process is simpler than other types of deep-inelastic collisions. The simplicity 
originates from the fact that fission is not sensitive to the fusion stage of a 
binary reaction whose understanding is far from clear. 

When mass and projectile energy increase, the possibility to isolate the 
ground-state to scission-point motion becomes problematic because of the in- 
creasing contribution of quasifission, a fission-like process occurring at angu- 
lar momenta J exceeding some value Jo where the fission barrier disappears. 
Quasifission in a clear cut sense has been observed in very heavy systems 
|T7| which do not have a fission barrier even at J — h. These experiments 
provided a motivation for the creation of the code HICOL which proved to 
be capable of reproducing the general features of the process ]18| . 

The Rayleigh function used in HICOL is built within the one-body dis- 
sipation model, proposed in Ref. It combines the wall formula in the 
mononuclear regime with an improved version of the so-called 'completed 



window formula' |20[ in the dinuclear stage. Calculations show that quasi- 



fission originates from the fact that the potential of the composite system 
for J pa J is very flat. Combined with strong friction this gives the system 
enough time to thermalize the relative velocity of the colliding nuclei and to 
relax the mass-asymmetry mode. 

In medium mass reaction systems, a quasifission configuration can exe- 
cute several rotations before scission, so that it is impossible to disentangle 
quasifission from fusion-fission by using fragment angular distributions as in 
heavier systems. Moreover, due to decreasing in angular momentum caused 
by the emission of few particles, a pocket can emerge in the initially flat 
potential, and a quasifission trajectory will go over into a fusion-fission one. 

The introduction of quasifission trajectories to fission-evaporation rou- 
tines is not an easy task. It requires answers on some principal questions. 
For example, in the fission case it is assumed 21] that, for all J within the 
fusion-fission angular momentum window, the system moves along the bot- 
tom of the fission valley of a non-rotating nucleus. Thus a question arises 
on whether the quasifission trajectories go along this path. Another related 
question is whether the evaporation model is applicable for particle emission 
along the quasifission trajectories. 
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The main purpose of the present work is to elucidate these questions. 
As a result we shall get a better perspective for a uniform description of 
fusion-fission and quasifission reactions accompanied by particle emission. To 
illustrate this opportunity, we test the compatibility of the charged particle 
clock with the neutron clock in calculations which effectively account for 
fusion-fission and quasifission trajectories. The study is performed in the 
A « 160, E x = 200-300 MeV region, where a quite complete set of data on 



multiplicities of fission-associated light particles is available [10,12 



The paper is organized as follows. In Sect. 2 we define the quasistatic path 
and in Sect. 3 we introduce the geometrical quantities needed to characterize 
the quasistatic and dynamic configurations of the system. These quanti- 
ties facilitate the comparison between the dynamic and quasistatic paths 
in Sect. 4. In Sect. 5 we perform Monte Carlo simulations of light-particle 
evaporation along dynamic trajectories and compare the results with avail- 
able experimental data. Our conclusions are drawn in the last section. 



0, J p{v)dv = A, (1) 



2 The Quasistatic Path 

The fundamental quantity in the construction of fission models is the total 
energy E{p} of the nucleus expressed as a functional of the spatial nucleon 
density p. The definition of this functional and its general properties are 
well documented (e.g., see Refs. p2|j23| , |24] , |25| ) . Given E{p] one can use the 
equations 

SE{p} 
bp 

where S/5p is the functional derivative, to find the ground state and the 
saddle point densities p gs and p s d- The difference between the corresponding 
energies, E{p sd } — E{p gs }, gives the fission barrier height. 

Since nuclei are leptodermous objects the terms 'density' and 'shape' are 
often used synonymously. As indicated first in Ref. PB"| , the so-called condi- 
tional equilibrium shapes may play an important role in nuclear dynamics. 
To define these shapes, one introduces a quantity Q characterizing the elon- 
gation of the nucleus. Q is a functional of the density and has the form 

Q{p} = J g(r)p(r)dr, (2) 

where q(r) is a known function. The density pq* of the conditional equilib- 
rium shape is the solution of equations (|l|) in the class of densities p restricted 
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by the condition Q{p} = Q* where Q* is a constant. As pointed out in Ref. 
f26fl , the sequence of conditional equilibrium shapes may acquire a physical 



meaning, for instance it can describe the process of fission if Q is changing 
adiabatically compared to all other degrees of freedom. It was conjectured 
in Ref. |26|] that the adiabaticity condition is likely to be satisfied in heavy 
systems where E{p} changes from the ground state to the saddle point by a 
few MeV only. 



Nowadays, there is a strong evidence |nj that the dynamical equation 
for fission of hot nuclei should take into account the coupling of Q with the 
single-particle degrees of freedom, which causes a slowing down of the fission 
time scale to values exceeding the single-particle times by a few orders of 
magnitude. This allows one to assume that at each instant of time the single- 
particle degrees of freedom are in thermal equilibrium under the constraint of 
the known Q and the respective velocity Q. Expressed in formal terms, this 
means that the entropy 5* of the system is maximal for the given constraints: 

— =0, Jp(r)dr = A, Q = Q\ Q = Q*. (3) 



Equations (|3p define a state of partial thermal equilibrium \T%\, or quasistatic 
state, for short. 

Denoting the intrinsic excitation energy by E x and using the Fermi gas 
formula for the entropy 



S = 2^a(E x -E{p} + E{p gs }), (4) 

one can find the density p of the quasistatic state by looking for the con- 
ditional minimum of E. This is justified by the fact that the variation of 
the level density parameter a with shape is very smooth in comparison to 
that of E p?J. By changing Q* one gets a sequence of quasistatic states 
which we call the quasistatic path. The fact that fission proceeds along the 
quasistatic path leads to significant simplifications in the formal description 
of this process. It also allows us to employ the statistical-model treatment 
of particle emission for the fissioning nucleus, because the single-particle de- 
grees of freedom are in thermal equilibrium at each point of the quasistatic 
path. 

Another process in which the quasistatic shapes can be useful is quasi- 
fission. The complete thermal equilibrium of the single-particle motion at a 
given shape may be reached during the reseparation stage of a quasifission 
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reaction if in the fusion stage the mass asymmetry mode and the relative 
velocity of the two colliding nuclei have relaxed. As a result, the quasifission 
trajectories in the reseparation stage will get on the quasistatic path. In 
practical terms, quasifission reactions are described by the one-body dissipa- 
tion model of heavy ion collisions |1| implemented in the code HICOL, while 
the conditional equilibrium densities can be calculated in the framework of 



the extended Thomas-Fermi (ETF) model of non-spherical nuclei I 27| , |28y29 |. 
With these two models we wish to verify whether the dynamic trajectories 
in quasifission reactions follow the quasistatic path. 

Our ETF calculations will be confined to mass and axially symmetric 
(about the z axis) prolate density distributions p(r, z) normalized to the mass 
number A, where r = y/x 2 + y 2 and x, y, z are the Cartesian coordinates. 
The energy density of the ETF model incorporates second-order gradient 
corrections with spin-orbit and effective mass terms which are very 
important in describing the nuclear surface. We have performed the ETF 



calculations using a realistic Skyrme interaction, namely SkM* ||24j| . From 
the ETF-SkM* functional, we obtain fully self-consistent nuclear densities 
by solving the associated variational Euler-Lagrange equations in cylindrical 



coordinates, imposing a given value of the quadrupole moment Q 2 p7| , p8| , |29| : 

roc rco r -, 

Q 2 = 2n / 2z 2 -r 2 p(r,z)rdrdz. (5) 



To account for nuclear rotation with angular momentum J we have included 
a rotational energy 

£rot = ^ (6) 

to the ETF energy functional. Here, / is the rigid-body moment of inertia 



7T 



OO /'OO 

2z 2 + r 2 







m p(r, z) r drdz , (7) 



where m is the nucleon mass. It was assumed that the spin axis is directed 
perpendicular to the symmetry axis of the compound nucleus. 

In HICOL the constant density approximation is used. Therefore, the 
nuclear density is completely determined by the profile function y(z) whose 
rotation about the symmetry axis z generates the nuclear surface. This 
poses difficulties in the comparison between dynamic and quasistatic den- 
sities, which can nevertheless be avoided by introducing some generalized 
characteristics of the nuclear densities. 
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3 Geometry of the Composite System 



For the sake of comparison between dynamic and quasistatic paths we in- 
troduce the elongation coordinate D mm and the neck coordinate -R nec k- In 
axially symmetric nuclei these quantities are given by 



roo poo 

A 



D mm = ~y I / z p(r, z) r drdz (8) 



Jo 
and 

2 r°° 

Rl ck = - rp(r,z = 0)dr. (9) 
Po Jo 

The elongation coordinate D mm defines the distance between the centers 
of mass of the two halves of the nucleus. It was used by Strutinsky |26| 



as a constraint operator in the integro-differential equation for the profile 
function y(z) of leptodermous nuclei. Our definition @ of the neck radius is 
obtained from the requirement that a nucleus with constant density po and 
a geometrical neck radius equal to -Rneck, has the same number of particles 
in the cross section of its neck as the nucleus having the distributed density 
p. In the following, po in Eq. (|9p will be identified with the one used in the 
code HICOL, namely p = A/(^-kRI) where Rq = 1.18 A 1 ^ 3 fm. From Eqs. 
(H) and (§) one finds that for a spherical nucleus with a constant density po, 
the elongation D mm is equal to |i?o an d that the neck radius is -R nec k = Ro- 
In HICOL the profile function y(z) of the composite system is parame- 
terized by two spheres smoothly connected by a second-order curve |30| . For 



m • 



these so-called Blocki shapes we define D mm as the distance between the cen- 
ters of mass of the two parts of the nucleus on both sides of a plane z = z- 
For z m we take the mean value of the left and right matching points (z m = 
for symmetric shapes). We identify -R noc k of the Blocki profile with y(z m ). 
For mass symmetric shapes, the Blocki profile function reads 

f Rj-{z + s/2f for - R x - s/2 < z < -z x , 
y 2 (z) = I a + pz 2 for - z x < z < z x , (10) 

[ R\-{z- s/2) 2 for Zl < z < s/2 + R x . 

Given the volume Vo of the nucleus, the parameters R\, zi, a and (3 can be 
expressed in terms of the two collective degrees of freedom s and a, where s 
is the distance between the centers of the spheres and 

Vo ~ 8vr^f/3 
a = Vo (H) 
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is a measure of the constriction of the system. 
Indeed, from Eq. flTT|) one obtains 



Ri — R 



a 



1/3 



;i2) 



where Rq is the radius of a spherical nucleus of volume Vq. If we equate the 
values of y(z) and its derivative on both sides of z = Z\ and require that the 
total volume of the shape generated by y(z) equals V , we get 



a 



Rl + 



z% - 



P = 7h 



s 

s 



with 



2-1 



G — — ^ 

4 



2 
2 1 



a 



s 6 1 



a 



so 



2R 1 



For the Blocki shapes defined by Eq. (|T0|), Eq. 
3 



yields 



A, 



s 
48 



1 _ G )3/2 ± S 2 R 2 + l SjR 3 + R A 

2 3 



(13) 
(14) 
(15) 

(16) 
(17) 



The value of R 2 eck coincides with y 2 (0) and according to (|i~0| ) is equal to a. 
Thus, using Eqs. flI3| ) and (|HS| ) we obtain 



neck 



R{ 



(18) 

allow one to 



Given R , Eqs. (0) and (pD together with Eqs. (p|) and ([ 
express £) mm and -R nec k in terms of s and a. 

In the following section, instead of the quadrupole moment, we use some- 
times the moment of inertia as a constraint operator. It should be noted that 
for mass symmetric Blocki shapes these quantities can be found analytically: 



'h(R 1 + h) 3 ( h\ h, 9 >9N „ 2 >9 ^ hzf 
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20 



(19) 
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nmp 



( J R 1 + /i) 3 /8 



R\ 




+ 3 [Kl )Zl+ 15 1+ 30 



(20) 



where /i = s/2. 



4 Dynamic and Quasistatic Paths 

We now describe calculations performed in order to compare the dynamic 
trajectories with the quasistatic path. The calculations are carried out in 
the A m 160 composite mass region, a region of continuous experimental 
efforts [p],p|,^|fl|7|fl|9| J10|J12l] . We start with the quasifission reaction following 



a 60 Ni+ 100 Mo collision at the beam energy E = 600 MeV. In Fig. 1 we dis- 
play the equidensity contour plots corresponding to the sequence of the ETF 
densities of conditional equilibrium for 160 Yb, which represents the compos- 
ite system in the collision. For each ETF density shown we have calculated 
the values of D mm and -R nec k by means of Eqs. (H) and (|9|). Inserting these 
D mm and -R nec k into Eqs. (|H]) and ( |T8"D and solving these equations with 
respect to s and a, we can prescribe the Blocki profiles to the density distri- 
butions. From Fig. 1 one can see that such profiles are close, in general, to 
the |po curves of the ETF density distributions. In dinuclear configurations, 
however, the nascent fragments predicted by the ETF model are somewhat 
flattened in the r direction compared to the spherical form. 

The densities shown in Fig. 1 were calculated accounting for the rotational 
energy with angular momentum J= 86 h. The reasons for this choice of J 
will become clear in the next section. The calculations indicate that the 
quasistatic path depends weakly on the angular momentum. When one goes 
from J=86 h to J=0 h the quasistatic D mm and -R nec k change at most by 
0.1 fm. In earlier ETF calculations p7| the energy E of the rotating nucleus 
along the quasistatic path was found to be very close to the sum of the 
energy of the non-rotating nucleus and the rotational energy computed with 
the moment of inertia of the latter. Moreover it was shown that the level 
density parameter is affected by rotation in a negligible amount if Qi is 
fixed. These findings were interpreted as an indication that rotation has a 
small influence on the nuclear densities calculated in conditional equilibrium. 
Our direct calculations of shape parameters along the quasistatic path agree 
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with this conclusion. 

We repeated the calculation of the J— 86 ft quasistatic path with the 
constraint on the moment of inertia /, instead of Q 2 - The path in the 
(.D mm , -Rncck) space turned out to be almost the same as the one we had 
found with the Q 2 constraint: the differences in -R n eck are no larger than 
~1%. Imposing the I constraint we not only obtain the same values for Q 2 
as with the Q 2 constraint, but also find that the hexadecapole moment Q4 
along the path is very similar. This means that the nuclear shapes must be 
equivalent with the / or Q 2 constraint, as well. 

The dynamic evolution of the shape of the 60 Ni+ 100 Mo composite system 
(at J=86 ft and E= 600 MeV) computed with HICOL is shown in Fig. 2. 
In units of 10 -21 s, the first stage (when neck fills in) takes about 0.2, the 
mononucleus lives about 30 and the scission stage lasts about 5. At the end 
of the first stage the individual temperatures and angular velocities of the 
two nuclei, predicted by HICOL, become approximately equal. The time 
dependence of D mm for different J values is depicted in Fig. 3, which shows 
a clear separation of the reactions into the fusion stage taking a fraction of 
10 -21 s and a much longer reseparation stage. The duration time of the latter 
strongly changes from one J to another. For example, this time decreases 
from 45 x 10~ 21 s to 15 x 10~ 21 s when J increases from 85 h to 100 h. 

Figure 4 shows on a (D mm , -R ncc k) plot how the dynamic trajectories with 
J=70, 75, 80, 85, 90, 95, 100, 105 ft are joining the quasistatic path deter- 
mined for J— 86 ft. The systems with J =70 and 75 ft terminate at the 
different points of this path (in the case of J=75 ft this happens after a slight 
rebound) while those with J =80-105 ft, having got on the quasistatic path 
shortly after rebound, proceed towards the scission point. In the vicinity of 
the scission point, they start to deviate progressively from the quasistatic 
path and yield a thinner neck for the same value of the elongation. We 
have found the dynamic paths, in the reseparation phase, to be rather stable 
against variations in the excitation energy and mass asymmetry. This follows 
from our dynamic calculations performed for the system 60 Ni+ 100 Mo at the 
beam energy of 1200 MeV (Fig. 5), and for 48 Ca on 112 Sn at the beam energy 
of 480 MeV (Fig. 6). 

Very small but regular deviations of the dynamic trajectories in their 
slowest phase from the quasistatic path are clearly observed in Figs. 4, 5 
and 6. They are probably related to the fact that the dynamical and varia- 
tional calculations involve different forces, and that in the ETF calculation 
we varied the whole density rather than just the nuclear profile. To verify 
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this assumption we calculated the quasistatic path for J= 86 H, with a Q2 
constraint, using the Yukawa-plus-exponential (YPE) forces which are 
employed in the code HICOL. We looked for the conditional minimum of 
the system energy in the class of mass symmetric Blocki shapes. The result 
is displayed on Fig. 6 with square symbols. According to this figure, the dy- 
namical paths in their slowest part practically coincide with the quasistatic 
YPE path. 

It is interesting to note that the YPE quasistatic path predicts the onset of 
the scission stage at noticeably smaller values of D mm than the Skyrme path. 
This is consistent with the fact that the saddle point configuration for the 
YPE force is more compact than for the SkM* force. The deformation energy 
for the YPE force reaches its maximum (of about 27.8 MeV) at D mm ~ 13.6 
fm, whereas the SkM* energy reaches its maximum (about 25.1 MeV) at 
-Dmm ~ 15.6 fm. These deviations manifest the scale of errors introduced 
into variational calculations by simple parameterizations of nuclear shapes. 



5 Particle Emission 

In the preceding section we have shown that the dynamic trajectories in 
the reseparation stage closely follow the quasistatic path. This means that 
statistical models of particle emission can be applied in this stage. Below, 
we describe the technique for light-particle evaporation calculations along the 
slow phase, which will allow us to perform comparisons with experimental 
data. 

Our procedure is based on a Monte Carlo simulation of particle decay 
chains in nuclei with a time-dependent shape. Given the excitation energy 
E x , angular momentum J, the dynamic path and the time t x for the beginning 
of the slow phase, we calculate the neutron (R n ), proton (R p ) and alpha parti- 
cle (R a ) emission rates. We assume that the emission times t e are distributed 
according to the exponential law exp(— R tot t c ), where R to t = R n + Rp + R a 
is the total emission rate. A specific value of t e is chosen using a genera- 
tor of exponentially distributed random numbers. Having sampled the type 
of the emitted particle in proportion to the weights R n / Rt ot , R p /Rt ot and 
Ra/Rtot, we find the average excitation energy u and root mean square an- 
gular momentum j rms = \J~p of the corresponding daughter nucleus and take 
the latter as the new decaying nucleus. Assuming its shape to be the same 
as that of the parent nucleus at t[ + t e , we simulate the next decay. The chain 
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of decays is terminated when the emission time exceeds the scission time of 
the nucleus and one proceeds to the next chain. 

The pre-scission multiplicities M V (J) are calculated as 



(21) 



where M v (i) is the number of particles of type v {y — n, p, a) in a decay chain 
classified with the index i. Here, M is the total number of decay chains for a 
given E x and J. The explicit expressions for R u , u, j 2 are summarized in the 
Appendix. They were obtained in the framework of the classical statistical 



531. As shown in 



model of particle emission from non-spherical nuclei 
Ref . |35j , the effects of the shape distortions on particle emission are treated 



by these formulas more rigorously than in heuristic models P<3|J57|J5^|J5P|,|^D 



The input parameters of these formulas are the effective separation en- 
ergies S^f calculated including deformation energies [11]], the level density 
parameter a, the height H, and the radius of the corresponding spherical 
barrier experienced by a particle. In the following calculations we employ the 
YPE values of S„ for v = n, p, a. In Fig. 7 the YPE values (solid lines) are 
compared with the ETF separation energies (dashed lines) on the quasistatic 
path in 160 Yb. The ratio I / Jo of the deformed nucleus moment of inertia to 
that of the spherical one is used as the coordinate along the path. 

The level density parameters a to be used later have been normalized at 
the spherical shape to the experimental value A/8. 8 MeV" 1 obtained in Ref. 
41f| . The shape dependence of a has been calculated with the YPE forces 
following the prescription of Toke and Swiatecki J42f. As seen from Fig. 8, 



these a are close, by magnitude and shape dependence, to the a values from 
the ETF method. The dependence of Vb and i?b for n, p, a on A and Z is 
parameterized in the same way as in Ref. 



43 



The energy of the emitted 
neutron entering the corresponding is replaced with its mean value ~ 2r, 
where r is the temperature of the daughter nucleus. 

Recently, Lou et al \ Y1\ measured multiplicities of light particles in fission 
reactions of 10 MeV/u 63 Cu+ 92 ' 100 Mo and 20 MeV/u 20 Ne+ 144 ' 148 ' 154 Sm. Our 
analysis of these data along with the data of Gonin et al |10] on 60 Ni+ 100 Mo 



is presented in Fig. 9 and Table 1. The measurements of multiplicities of n, 
p, a in the latter work were performed at 9.2 and 10.9 MeV/u. The data we 
analyze are obtained by interpolation of the reported values to 10 MeV/u. 
In the analysis we took into account that in these reaction systems an 
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appreciable amount of particles escape from the system during the pre- 
equilibrium stage. The measurements of linear-momentum transfer from pro- 
jectile to target allowed Lou et al ]I2] to estimate that the mass and charge 
removed in this stage are (6A,5Z)= (8,4), (8,4), (6,3), (7,3) and (6,3) in the 
reactions 63 Cu+ 92,100 Mo and 20 Ne+ 144 ' 148 ' 154 Sm, respectively. Since experi- 
mental information on 5 A and SZ for the 60 Ni+ 100 Mo system is not available, 
we used the values (SA, 5Z) = (8,4) in the closest system, 63 Cu+ 100 Mo. Lou 
et al [12| estimated the initial excitation energies of the equilibrated com- 



pound nuclei to be E x = 227, 267, 289, 277 and 282 MeV in the reactions 
of 10 MeV/u 63 Cu+ 92 ' 100 Mo and 20 MeV/u 20 Ne+ 144 - 148 ' 154 Sm, respectively. 
For the system of 10 MeV/u 60 Ni+ 100 Mo, the linear interpolation between 
the excitation energies E x = 251 MeV and 293 MeV at 9.2 MeV/u and 10.9 
MeV/u, respectively, obtained in Ref. [[Hj, results in E x = 271 MeV. 



This information on 5 A, 5Z and E x was used as input in our pre-scission 
(equilibrium) multiplicity calculations. For all systems we used t\ = 1.1 x 
10 _21 s and simulated M = 200 chains of decays. For each system we took 
only one HICOL trajectory, namely the trajectory whose contact-to-scission 
time t cs is closest to 40 x 10~ 21 s. This condition results in J= 78, 86, 86, 
85, 88, 93 h in the reactions of 63 Cu+ 92 Mo, 60 Ni+ 100 Mo, 63 Cu+ 100 Mo at 10 
MeV/u and 20 Ne+ i44,i48,i54g m at 2 q MeV/u, respectively. The slow stage 
of the so-chosen ('effective') trajectories lasts some 35 x 10~ 21 s. This is 
consistent with the fission time scale of (35 ± 15) x 10~ 21 s deduced from a 
systematic study of pre-scission neutron multiplicities in 27 fission reactions 
induced by 16 ' 18 0, 40 Ar and 64 Ni on targets with A = 141-238 0. 

In the post-scission emission calculations, the thermal energy of the com- 
posite system at the moment of scission was shared between the comple- 
mentary fragments in proportion to their masses. The spins of the fragments 
(about 6 h per fragment) were taken from HICOL output. To find A and Z of 
the primary fragments, we used the calculated n, p and a pre-scission multi- 
plicities (see Table 1). The calculated post-scission multiplicities of neutrons 
in all reactions but 63 Cu+ 92 Mo are confined between the value of 3.7±0.4 
measured in the system of 16 0+ 154 Sm at E x = 206 MeV || and the value of 
4±1.1 for 60 Ni+ 100 Mo at 10 MeV/u which follows from interpolation of the 
9.2 and 10.9 MeV/u data (3.6 ± 1 and 4.5 ± 1.2, respectively) of Ref. 0. 

With accounting for post-scission emission, which is essential, in fact, 
only for neutrons, the total (equilibrium) multiplicities appear to be within 
the likely systematic uncertainties of the experimental points. The only 
noticeable exceptions are the total proton multiplicities in 63 Cu+ 92 Mo and 
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60 Ni+ 100 Mo, when the measured values are smaller than the calculated ones 
by factors of 1.6 and 2.1, respectively. With the exception of these two data 
points, the overall agreement of the calculations with the rest of the data 
signifies a consistency of the neutron with the charged particle clock. The 
low proton multiplicities in 60 Ni+ 100 Mo reactions observed in Ref. fl(| have 



motivated further experimental studies. In a recent work, Charity et al |44 
studied the nearby system 64 Ni+ 100 Mo at a similar excitation energy and 
found 2-3 times greater multiplicities of /kszon-associated p and a compared 



to those of Ref. 10 



It should be noted that no fitting parameters were used in our analysis. 
The employment of the transmission coefficients for p and a from Ref. [45| 



would destroy the quality of the description by strongly enhancing the a 
emission. For instance, the pre-scission multiplicities in the 60 Ni+ 100 Mo sys- 
tem become 6.73, 2.29 and 1.81 for n, p and a, respectively, instead of 7.27, 
2.11 and 0.98 obtained with the transmission coefficients from Ref. ||43||. 



6 Discussion and Conclusion 

To shed more light on the role of quasifission trajectories in the reactions of 
our study, it is useful to estimate the fusion-fission J-window and compare 
it with available experimental information on the angular momenta, associ- 
ated with the evaporation residue and fission cross sections. Nuclei emerg- 
ing at the end of the evaporation cascades undergo fission if their angular 
momentum is confined between the angular momentum where the fission 
barrier height -Bf(J) equals the neutron separation energy and the angular 
momentum where -Bf(J) vanishes 0,0. To reconstruct the corresponding 



J-window in the beginning of the evaporation cascades, one has to account 
for the angular momentum removed by pre-scission particles. 

In the case of 60 Ni+ 100 Mo at 10 MeV/u, we find that the fusion-fission 
J-window at the end of the evaporation cascades is 57-82 h, on the average. 
Light particles evaporated along the effective trajectory, remove on the aver- 
age about 21 h. Therefore, we estimate the corresponding J-window in the 
beginning of the equilibrium emission stage, i.e. for 152 Dy, to be 78-103 h. 
Since B{{ J) for 152 Dy vanishes at 83 h, trajectories within this window, with 
the exception of those with J= 78-83 h, belong to quasifission. 

It is interesting to note that the width of the so-defined J-window is close 
to the 'total' (fusion-fission plus quasifission) J-window widths in the closest 
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systems to ours, where data are available. A fission J-window of 70-103 
% and 49-67 % is implied in studies of 40 Ar+ 109 Ag at 8.4 MeV/u [g§] and 
20 Ne+ 159 Tb at 16 MeV/u ^9[, respectively. Thus we may conclude that 
the majority of fission-evaporation reactions start on quasifission trajectories 
which after losing angular momentum end up on trajectories in the potential 
with a non-zero fission barrier. 

To summarize, the present work was motivated by the desire to estimate 
the perspective for the inclusion of quasifission into fission-evaporation codes. 
Towards this aim we tested whether the quasifission trajectories, in the re- 
separation stage, follow the fission path. Our calculations were performed in 
the A ph 160 region at the bombarding energy of 10-20 MeV/u. 

For the sake of comparison between different shapes (assumed to be axi- 
ally symmetric) we use the two-dimensional space of the collective variables 
describing elongation and constriction. The elongation is characterized by the 
distance between the two halves of the nucleus. The constriction is described 
by the neck radius defined in a way applicable for distributed densities. In 
this collective space, we first check the sensitivity of the quasistatic path to 
the angular momentum of the system and to the form of the constraint oper- 
ator. The quasistatic paths found for different angular momenta (including 
those exceeding the critical angular momentum for fission) are practically 
indistinguishable from the one with J= H. The moment of inertia as a con- 
straint operator was found to generate a sequence of shapes which coincides 
with the one obtained using the constraint on the quadrupole moment. 

The same space of collective variables was used to compare the sequences 
of shapes predicted by the HICOL code at different values of the entrance- 
channel angular momentum. The parts of the quasifission trajectories corre- 
sponding to the slow stages of the evolution show the eventual convergence to 
the quasifission trajectory with the minimal J. This latter is found to follow 
closely the quasistatic path obtained with ETF model. This close coincidence 
of dynamic trajectories with the quasistatic path was found to occur in a wide 
range of bombarding energies and for quite different entrance-channel mass 
asymmetries. 

Since the dynamics of quasifission reactions is well described by qua- 
sistatic paths, at least during the slowest phase, the statistical evaporation 
model can be applied for the description of particle emission from such sys- 
tems. This follows from the fact that at each point of the quasistatic path 
the system reaches thermal equilibrium. 

Using single effective dynamical trajectories whose slow stage lasts some 
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35 x 10~ 21 s, it was made possible to reproduce reasonably well experimental 
data on pre- and post-scission multiplicities of neutrons and total multiplici- 
ties of protons and a-particles emitted from thermally equilibrated systems. 
This agreement was achieved without any ad hoc statistical model parameter 
adjustments and shows a consistency between the neutron and charged par- 
ticle clock. The duration time of the slow stage of the employed dynamical 
trajectories was found consistent with the results of systematic studies. 
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Appendix 

In the following we outline the formalism we used for the calculation of decay 
rates for statistical particle emission from equilibrated compound nuclei. The 
excitation energy, angular momentum, deformation energy and moment of 
inertia of the parent nucleus are denoted as E x , J, E^ei and I x , respectively. 
Given these quantities, the thermal energy of the parent nucleus is defined 



by 



— Ex — 777 ^def- 



x 



Its reduced level density is 



1 



r exp [2y/a x q x \ , 



where 
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and a x is the level density parameter. Similar formulas are used for the 
reduced level density w(q) of the daughter nucleus with thermal energy q. 
Its moment of inertia and level density parameter are denoted by / and a. 

The mass, spin and effective separation energy of the emitted particle 
are denoted as m„, s u and Sf 1 , respectively. Let -Zmatter be half a length of 
matter distribution in the deformed shape, and let R^ and -R mat tcr be the 
barrier radius and matter radius for the spherical shape. Then half a length 
of the figure generated by the barrier line is postulated to be 

ZO — ^matter + (Rb ~ R matter)- 

Given the matter profile function y = y(z) and the Coulomb potential &{z) 
along this profile, the barrier line p(z) and the potential barrier U (z) along 
this line were calculated from 



= Ky (±) , m = |* (±) , 



where K = zo/z mat ter is the scaling factor, Vt, is the s-wave potential barrier 
in the spherical nucleus and $o is the Coulomb potential on the edge of the 
spherical matter distribution. 

The key characteristics of the residual nucleus are its average thermal 
energy q and temperature r corresponding to this thermal energy. These 
quantities are calculated from 



E„ — S 



eff 



E, 



dcf 



J 2 

21 - 



a V \aj a 



With these definitions the particle emission rate is given by 



2s u + 1 22 u(q) 
R v = — m u z T — - ( — ^ exp 



2vr 



\\E x -Sf-E Aci - 3 --q 



G(r,b), 



where 



G(T,f>) = j^d(yf + ^(^l\ exp 



T ^2 
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m v zl (J\ 2 p((z ) 



72 _ 72 :_,/ /lll6 ^ 7 - / ' ) 



2r V// ' ' v y z 

and Iq(x) is a Bessel function of the first kind of imaginary argument. 

The averaged square of the angular momentum of the daughter nucleus 

is 

Given this quantity and q, the average excitation energy of the daughter 
nucleus reads _ 

f 

u = q + — + M dc f, 

where i^ef is the daughter nucleus deformation energy. 

The above expressions give approximate values of R v , j 2 and u because 
they use approximate values of q and r. To correct these latter quantities we 
find a new value of q using the formula 

H del 2/ ^ 

at the initial value of r and insert it into the formula for r. This procedure 
was iterated until convergence of r was attained. 
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Table 1. Comparison between experimental and calculated light-particle 
multiplicities for the reaction systems studied in Refs. |I| and fL2| . The 
neutron number iV of the composite system after the pre-equilibrium emission 
stage is indicated. From left to right the systems are 10 MeV/u 63 Cu+ 92 Mo, 
60 Ni+ 100 Mo, 63 Cu+ 100 Mo and 20 MeV/u 20 Ne+ i44,i48,i54g m 



N 


80 


86 


88 


89 


92 


99 


n pre 


3.75 


7.27 


7.46 


7.65 


8.2 


10.87 


^post 


1.79 


3.84 


3.87 


3.82 


3.89 


3.97 


j^tot 


5.54 


11.11 


11.33 


11.47 


12.09 


14.84 


j^tot,exp 


6(1) 


11.25(2.7) 


10.9(1.5) 


10.3(1.3) 


12.2(1.6) 


16.9(1.2) 


p prc 


2.97 


2.11 


1.97 


2.65 


1.85 


0.82 


ppost 


0.18 


0.06 


0.05 


0.07 


0.05 


0.02 


ptot 


3.15 


2.17 


2.02 


2.72 


1.9 


0.84 


ptot,exp 


2.02(0.3) 


1.03(0.15) 


1.6(0.3) 


2.14(0.3) 


1.83(0.5) 


0.97(0.5) 


a pre 


1.17 


0.98 


0.93 


1.24 


1. 


0.66 


Q/POSt 


0.03 


0.02 


0.01 


0.02 


0.01 


0.01 


a tot 


1.2 


1 


0.94 


1.26 


1.01 


0.66 


Q,tot,exp 


1(0.3) 


0.77(0.11) 


1.69(0.4) 


1.3(0.3) 


0.7(0.4) 


0.63(0.4) 
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Figure Captions 



Fig. I. Equidensity contours for 160 Yb calculated by the ETF method at 
the indicated values of the quadrupole moment, with angular momentum J— 
86 h. From outside to inside the lines represent contours of constant density 
p = 0.1, 0.3, 0.5, 0.7, 0.9 and 1.1, in units of pq. The dashed curves represent 
the Blocki profile with the parameters s and a calculated in terms of -D mm 
and .R n eck obtained from the corresponding ETF density distribution. 

Fig. 2. Time evolution of the 60 Ni+ 100 Mo system from the code HICOL 
at E = 600 MeV and J= 86 H. The time is indicated in units of 10~ 21 s. 

Fig. 3. Time dependence of D mm for the reaction 60 Ni+ 100 Mo at E = 600 
MeV and J= 79-105 h predicted by the HICOL code. 

Fig. 4. Dynamic (for J = 70, 75, 80, 85, 90, 95, 100, 105 K) and quasistatic 
(for J= 86 K) paths in the 60 Ni+ 100 Mo collision at 600 MeV are shown as 
solid and dashed lines, respectively. The spherical shape is shown as a cross. 
Numbers along the curves indicate the time in units of 10~ 21 s for the dynamic 
path at J=86 Ti. The insert shows the final stage of the fusion trajectories 
with J=70 h and J=75 h. 

Fig. 5. Dynamic paths (for J = 70, 75, 80, 85, 90, 95, 100, 105 K) in the 
60 Ni+ 100 Mo collision at 1200 MeV are shown as solid lines. The quasistatic 
path for J= 86 Ti is represented by the dashed line. The spherical shape is 
shown as a cross. 

Fig. 6. Dynamic paths (for J= 70, 75, 80, 85, 90, 95, 100, 105 h) for 
the collision 48 Ca+ 112 Sn at 480 MeV are shown as solid lines. The dashed 
line represents the quasistatic path corresponding to the self-consistent ETF 
variational calculation with the SkM* force. The closed square symbols show 
the quasistatic path found for the YPE forces in the space of Blocki shapes. 

Fig. 7. The effective separation energies of neutrons, protons and alpha 
particles along the quasistatic path in 160 Yb. The solid lines show the cal- 
culation for the dynamic shapes with YPE forces. The dashed lines indicate 
the ETF effective separation energies. 

Fig. 8. The level density parameter along the quasistatic path. The solid 
line represents the a values normalized to A/8.8 MeV -1 for the spherical 
shape with the Toke-Swiatecki shape-dependent factor based on YPE forces. 
The dashed line indicates the ETF calculation. 

Fig. 9. Equilibrium total multiplicities of neutrons, protons and alpha 
particles (pre-scission plus post-scission) as a function of the neutron num- 
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ber iVof the emitting system (after the pre-equilibrium stage). Experimental 
data points from Refs. |L2] and |L0| are shown by filled squares and crosses, 
respectively. The short-dashed lines represent the calculated pre-scission 
(equilibrium) multiplicities. The solid lines show the total calculated mul- 
tiplicities. In the top (neutron) panel, the open circles represent the total 
measured multiplicities accounting for pre-equilibrium emission [12 . From 
left to right the systems are 10 MeV/u 63 Cu+ 92 Mo, 60 Ni+ 100 Mo, 63 Cu+ 100 Mo 
and 20 MeV/u 20 Ne+ 144 ' 148 ' 154 Sm. 
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Figure 2 
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